Low lying states in ^B 
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Excitation functions of elastic and inelastic '^Be+p scattering were measured in the energy range 
between 1.6 and 2.8 MeV in the cm. An R-matrix analysis of the excitation functions provides 
strong evidence for new positive parity states in *B. A new 2^ state at an excitation energy of 2.55 
MeV was observed and a new 0"*" state at 1.9 MeV is tentatively suggested. The R-matrix and Time 
Dependent Continuum Shell Model were used in the analysis of the excitation functions. The new 
results are compared to the calculations of contemporary theoretical models. 
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Light nuclei are of great importance in modern nu- 
clear physics as their structure provides a link between 
nucleon-nucleon interactions and macroscopic nuclear 
many-body dynamics. Generally, properties of stable 
light nuclei, including level schemes, are reproduced 
rather well by the so-called ab initio methods that start 
from the basic interactions of nucleons jl, 2]. However, 
the neutron deficient isotope of boron, ^B, and its mir- 
ror nucleus ®Li, provide an interesting exception. Most 
ab initio calculations, including one from more than 10 
years ago 3 1 predict more positive parity states below 4 
MeV than what has been observed experimentally. These 
"missing" states are proposed to have a relatively simple 
structure, with large spacing between the levels and at 
low excitation energies. Consequently, it is rather sur- 
prising that these levels have not been observed to date. 
The importance of the ^Be(p,7)^B reaction for under- 
standing the solar neutrino flux is another stimulus for 
taking a closer look at the ^B structure. Previously un- 
accounted for low-lying states in *B may alter the the- 
oretical extrapolation of the ^Be(p,7)*B S-factor. The 
main objective of this work is an experimental search for 
these proposed low-lying levels in *B. 

The level structure of ^B and ^Li below 4 MeV is shown 
in Figure[TJ The first and second excited states, the 1+ at 
0.7695 MeV and the 3+ at 2.32 MeV, have been observed 
in numerous experiments, reviewed in f4| . The broad neg- 
ative parity state at ^-^3.0 MeV was first suggested in Ref. 
[3| and later identified as a 2^ state at 3.5±0.5 MeV in 
Ref. [6j . An excitation function of ^Be-|-p resonance elas- 
tic scattering was measured in both of these works. The 



most recent Be+p measurement, performed by H. Yam- 
aguchi, et al., [7], confirmed the 2~ state and determined 
the excitation energy and width of this state with better 



precision, 
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FIG. 1. (Color online) Experimental level scheme of *B and 
*Li at excitation energies below 4 MeV. New states are shown 
in red. The dash-dotted line indicates that the state is tenta- 
tive. 



positive parity states in the ^Be-|-p elastic scattering ex- 
citation function below 4 MeV have been observed in any 
of these studies. However, it is possible that the missing 
states still exist but contribute very little to the elastic 
excitation function due to the strong decay branch to the 
1/2" first excited state of "^Be at 0.43 MeV. This possibil- 
ity was considered by D. Halderson in the framework of 
the Recoil Corrected Continuum Shell Model (RCCSM) 
^] , and the suggestion was made to use inelastic ^Bc+p 
scattering to search for the missing states. It is interest- 
ing to point out that some evidence for the 2+ state at 
3.0 MeV was also presented in Q based on the RCCSM 
analysis of the ^Be-t-p elastic scattering excitation func- 
tion measured in [6] . In the present work, the excitation 
functions of the ^Be+p elastic and inelastic scattering 
were measured simultaneously and a consistent R-matrix 
analysis of both excitation functions was performed. 

The nucleus ^B is weakly bound with a proton sepa- 
ration energy of only 137 keV. All of the excited states 
of this nucleus are in the continuum. The recently devel- 
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FIG. 2. (Color online) Scatter plot of kinematic coincidence 
between detected protons and ^Be ions. Regions which cor- 
respond to elastic and inelastic scattering are labeled. The 
schematic view of the experimental setup with a sandwich of 
three Micron Semiconductor S2 detectors for the light recoils 
and a single S2 detector for the heavy recoils is shown in the 
inset. 



oped Time Dependent Continuum Shell Model (TDCSM) 
approach Q bridges the reaction-structure gap. Within 
this approach the cross sections for elastic and inelastic 
nucleon scattering can be calculated directly from the 
nuclear effective Hamiltonian. Influence of the contin- 
uum on the wavefunctions of the populated resonances 
is treated self-consistently and the number of free pa- 
rameters is greatly reduced. Once residual interactions 
are chosen, only one free parameter (excitation energy) 
remains for each resonance (compared to five in two- 
channel R- matrix approach). TDCSM analysis of the 
resonance scattering data is presented here. It was used 
not only as a stand-alone tool but also as a logical start- 
ing point for the subsequent R-matrix fit. 

The experiment was carried out at the John D. Fox Su- 
perconducting Accelerator Laboratory at Florida State 
University. A radioactive beam of ^Be was produced us- 
ing the ^II(^Li7Be)n reaction. The primary ^Li beam 
was accelerated by a 9MV SuperFN Tandem Van de 
Graaff accelerator followed by a LIN AC booster. The 
primary target was a 4 cm long hydrogen gas cell with 
2.5 /Ltm Havar entrance and exit windows. The gas cell 
was cooled by liquid nitrogen and had a gas pressure of 
390 mBar. The radioactive nuclear beam facility RESO- 
LUT was used to separate ^Be from other reaction prod- 
ucts and the primary beam. Two ^Be beam energies were 
used in this experiment: 22.0, and 18.5 MeV. The typical 
intensity of the ^Be beam was 10^ pps. The composition 
of the beam was 70% ^Be with 30% ^Li contaminant. 

A sketch of the experimental setup is shown in the in- 
set of Figure [51 A solid polyethylene (C2H4) target of 
thickness optimized for the given beam energy (see de- 
scription below) was used. A set of three annular Micron 
Semiconductor sihcon strip detectors (S2 design) for the 
proton recoils were positioned 5, 6 and 7 cm downstream 
from the target. Another S2 detector for the ^Be recoils 
was positioned 24.5 cm from the target. The S2 detector 



has annular geometry and consists of 16 segments and 
a side of rings that allow for the scattering angle of the 
products to be determined. The first in the set of three 
proton detectors was a AE detector with thickness 65 
fim. It was only used in the initial part of the experi- 
ment to verify that correct identification of light recoils 
can be achieved. The other two proton detectors and the 
^Be detector were 500 ^m each in thickness. 

The target thickness was optimized for maximum en- 
ergy losses of the ^Be ions in the target while ensuring 
that all ^Be recoils make it out of the target with enough 
kinetic energy left to be detected in the downstream S2 
detector. Kinematic coincidence between protons in the 
array of three S2 detectors and the ^Be recoils in the 
downstream S2 detector was then used to identify the 
scattering events. The time between the events in the 
proton and ^Be detectors was measured with time res- 
olution of about 3 ns to eliminate random coincidence 
background. Elastic and inelastic scattering processes 
can be distinguished because the complete kinematics of 
the event is measured. More detailed description of the 
experimental technique can be found in [10| . 

The 2D scatter plot for the kinematic coincidence be- 
tween protons and ''Be is shown in Figure [H The kine- 
matic loci which correspond to elastic and inelastic scat- 
tering processes are labeled and outlined with contours. 
Polyethylene target thicknesses used in this experiment 
were 2.5 and 1.5 mg/cm^ for the 22 and 18.5 MeV beam 
energies respectively. In addition, a separate run at 18.5 
MeV of ^Be beam energy was performed with a slightly 
thicker (2 mg/cm^) target to extend the measured ex- 
citation function to lower energies without changing the 
energy of the beam. This use of a thicker target comes 
at the price of losing coincidences between the highest 
energy protons and the ^Be recoils because they are pro- 
duced at the very beginning of the target and never make 
it through. Only the lower energy part of this spectrum 
was used in the analysis. 

Figure [3] shows excitation functions of elastic (top 
panel) and inelastic (bottom panel) scattering of ^Be-|-p 
measured in three different runs. The open circles corre- 
spond to the run at 18.5 MeV of ^Be with a 1.5 mg/cm^ 
target, the solid circles are 18.5 MeV of ^Be with a 2 
mg/cm^ target data and the stars are 22 MeV of ''Be 
with 2.5 mg/cm^ target data. Binning of 4° in the lab 
frame was used. Absolute normalization of the cross sec- 
tion was done using the known excitation functions for 
the ''Li-l-p elastic scattering [11]. These excitation func- 
tions were extracted from the experimental data using 
the same procedure as for the ''Be-l-p elastic scattering. 
Excitation functions extracted from our data agree well 
with the differential cross section for the elastic and in- 
elastic scattering of ''Be-l-p measured at several energies 
of ''Be using the conventional thin target approach and 
reported by U. Greife, et al., J12i |. 

A striking feature of our data is that the cross section 
for inelastic scattering is very large (^30 mb/sr at an 
excitation energy of 2.5 MeV). Two channel, multi- level 
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FIG. 3. (Color online) Excitation functions of Be+p elas- 
tic (top) and inelastic (bottom) scattering. Data from three 
different runs are shown. The open and solid circles corre- 
spond to the run at 18.5 MeV of ^Be with a 1.5 mg/cm^ 
and 2 mg/cm^ targets respectively, and the stars are at 22 
MeV with a 2.5 mg/cm^ target. The black solid curve is 
an R-matrix fit with only the known 3^ and 2~ states at 
2.3 and 3.5 MeV respectively. The red short-dashed line in- 
cludes contribution of the higher lying 1^ states assumed at 
3.0 MeV. The blue dash-dotted line shows the 1^ state shifted 
to 2.5 MeV and the green long-dashed line also includes the 
1~ state introduced at 5 MeV. The black dash-double-dotted 
curve in the bottom panel shows the excitation function of 
^Be(p,p')^Be(0.43) reaction from the direct mechanism. 



R-matrix analysis clearly indicates that it is not possible 
to explain this high inelastic cross section if only known 
states in ^B are considered (Figure [3]). This failure can 
be understood from the following simple considerations. 
The first excited 1"*" state at 0.77 MeV is too narrow 
(35.6 keV) to have any significant impact on the exci- 
tation functions at energies above 1.5 MeV. The second 
excited state, S"*" at 2.32 MeV, can only decay to the 
1/2" first excited state of ^Be with orbital angular mo- 
mentum £ = 3. Therefore, even if the corresponding re- 
duced width is large, the inelastic partial proton width, 
Tpi — 2Pi{kR)'-f'^ , would still be small compared to the 
elastic partial proton width due to a small penetration 
factor for high angular momentum decay. Hence, the 
cross section for population of the first excited state in 
^Be due to the 3+ resonance in *B, determined by the 
^p^p' /^ tot ratio, is small. The same is true for the broad 
2~ state in ®B at 3.2 MeV as it can only decay to the 
first excited state in ^Be with angular momentum £ = 2 
while decay to the g.s. proceeds with £ = 0. The black 
solid curve in Figure [3] shows the results of an R-matrix 
calculation with only the previously known 1"*", 3"*" and 
2~ states at 0.77, 2.32 and 3.7 MeV with reduced widths 
parameters evaluated using the TDCSM approach and 
known total widths of these states. (Excitation energy 
and width of the 2~ were adjusted slightly to produce 
a better fit.) It is clear that while the elastic scattering 



data are well reproduced, the inelastic scattering data 
cannot be explained by the known states. 

An attempt has been made to reproduce the observed 
p-|-^Be inelastic scattering excitation function without in- 
troducing new resonances in ^B but assuming a direct 
excitation of the ^Be first excited state in p-f-^'Be scatter- 
ing. In this case the reaction does not proceed through 
the population of resonances in ^B and cannot be eval- 
uated using R-matrix approach. The calculations were 
performed using the coupled-channels approach with the 
computer code FRESCO 0. The black dash-double- 
dotted curve in the bottom panel of Figure |3] shows the 
excitation function of the '^Be(p,p')^Be*(0.43) inelastic 
scattering at 146° due to the direct mechanism. It is 
clear that unless additional resonance(s) are introduced 
direct excitation cannot be responsible for the large in- 
elastic scattering cross section observed experimentally. 
Details of this calculation will be published elsewhere. 

Based on the level scheme of ^Li it is natural to try 
to introduce the second 1+ state in *B at an excita- 
tion energy around 3 MeV. Reduced widths for this state 
were chosen according to the TDCSM calculations car- 
ried out with Cohen-Kurath interaction (CKI) [lj|. It 
was verified that these reduced widths reproduce the 
known width of this state in ^Li (~1 MeV). The TD- 
CSM predicts that this state has a substantial inelastic 
partial width. The red short-dashed curve in Figure [3] 
shows the effect of the 1+ state at 3.4 MeV on the fit. 
While the elastic excitation function is fitted perfectly, 
the inelastic cross section is still underestimated, even if 
this state is shifted to 2.5 MeV, where inelastic scattering 
has its maximum cross section (blue dash-dotted curve in 
Figure [3]) . Finally, in an attempt to increase the inelas- 
tic cross section without using new states below 3 MeV 
we introduced a "background" state, the 1~ at '--^5 MeV. 
This state is a spin-orbit partner to the known 2" state 
and splitting between these states should be 1-2 MeV 
based on TDCSM calculations. This state can decay to 
the first excited state of ^Be with £ = 0, therefore it may 
contribute significantly to the inelastic cross section. The 
reduced widths for the 1~ state were evaluated using the 
TDCSM (WBP [ill residual interactions were used). As 
expected, the 1~ state increases the inelastic cross sec- 
tion overall (green long-dashed curve in Figure |3]). But 
even with this state included the inelastic cross section 
cannot be reproduced. 

The "extra" low- lying states in *B predicted by ab 
initio and shell model calculations have spin-parity as- 
signments 0+, 1^ and 2+. Infiuence of the 1+ state has 
already been discussed. It was found that introducing 
a new 2+ state placed at 2.55 MeV, reproduces both 
the magnitude and angular dependence of the observed 
peak in the inelastic cross section while keeping the elas- 
tic excitation function in agreement with the experimen- 
tal data. This is shown as green dash-dotted curve in 
Figure m where ^Be-|-p elastic (panels a-c) and inelastic 
(panels d-f ) scattering excitation functions measured at 
three cm. angles are shown. However, even with this 
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FIG. 4. (Color online) Final R-matrix fit of the experimental data, with the new 2^ and 0^ states. ^Be+p elastic scattering 
data at three cm. angles are shown in the top panels (a-c) and inelastic data are shown in the bottom panels (d-f). Solid blue 
line is the best fit, dashed red line is a fit with 1^ at 1.9 MeV instead of 0^ and green dash-dotted line is a fit without the 0^. 
Dashed-double-dotted cyan line in the bottom panel is TDCSM calculations with the known states in B and the new 0^ and 
2+ at 1.9 and 2.55 MeV respectively. 



TABLE I. Parameters of resonances in *B from the R-matrix 
best fit. States in parenthesis are unknown states outside of 
the measured excitation energy range. They provide essen- 
tial "background" through low energy tails. New states in *B 
introduced in this work indicated with superscript (a). Un- 
certainties correspond to la. 
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new 2+ state the cross section for inelastic scattering be- 
low 2.3 MeV is still too low. The 2+ state should have a 
relatively small width (360 keV) to fit the observed peak- 
like structure in the inelastic excitation function at 2.5 



MeV, and its influence below 2.3 MeV is small. The only 
predicted state which has not been considered is the 0+. 
Introducing a new 0+ state at the excitation energy of 1.9 
MeV with a width of 610 keV allows the inelastic scat- 
tering data to be fit below 2.3 MeV without destroying 
the fit to the elastic scattering data (solid line in Figure 
21). It was verified that a 1"*" spin-parity assignment for 
this new state is not possible as it will ruin agreement 
with the elastic scattering data (red dashed line in Fig- 
ureS]). An important distinction between the O"*" and the 
2+ states has to be made. While existence of the 2+ state 
is hard to dismiss, the case for the 0+ state is somewhat 
weaker and further investigation is warranted. 

An important role in guiding our analysis and pro- 
viding support for the 0"*" state was played by the 
TDCSM calculations Q. The excitation function for 
^Be(p,p')^Be(0.43) inelastic scattering calculated using 
the TDCSM approach is shown in the bottom panels of 
Figure S] as a dash-double-dotted cyan curve. The only 
free parameters in these calculations are the excitation 
energies of the states. All known states in ^B and the 
new O"*" and 2+ states at 1.9 and 2.55 MeV were taken 
into account in these calculations. The CKI residual in- 
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FIG. 5. (Color online) Comparison of the experimental data 
to the predictions of the ab initio models. The subset "Exp" 
shows experimental data. The new states in *B observed in 
this experiment are shown in red. The dashed line is a 1"*" 
state inferred from level structure of ^Li based on mirror sym- 
metry. The subset "NCSMl" is the result of the most recent 
NCSM calculations with the CD-Bonn 2000 potential [3. 
The subset "GFMCl" is the earher GFMC result [l3| with 
AV18/UIX Hamiltonians (values are for *Li and excitation 
energies for the last two states, shown as green dash-dotted 
lines, were calculated using Variational Monte Carlo rather 
than GFMC method). The subset "GFMC2" is the more re- 
cent GFMC prediction [1] with AV18/IL2 Hamiltonians (val- 
ues are for **Li). The subset "NCSM2" is the earlier NCSM 
result j3]- 



teraction [ij] was used and the corresponding states were 
shifted to their experimental locations. Note that the 
cross section at ^2 MeV is well reproduced by the 0+ 
state. This can be considered as an additional argument 
in favor of the 0+ state at 1.9 MeV in *B. At ~2.5 MeV 
the TDCSM cross section is determined by the 2+ state 
and it is lower than observed experimentally. This is 
due partially to the absence of the negative parity states 
(specifically the 1~ state) in the applied TDCSM model 
space. Details of the TDCSM analysis will be published 
elsewhere. 

Table U shows the best fit parameters of the 7 states, 
which were introduced to describe both elastic and in- 
elastic data simultaneously. Three of these states, the 
first 1+, 3+ and 2~ were known. The two states shown 
in parenthesis are unknown "background" states. The 
initial fit parameters for these states were taken from the 
TDCSM calculations and also (in the case of the 1+ at 
3.4 MeV) from the level structure of the mirror nucleus, 
^Li. The new low-lying O"*" and 2+ states were introduced 
to fit the measured cross section for p-l-'^Be inelastic scat- 
tering. The location of these states in ^Li is not known. 



calling for a new experimental effort to locate them in 
8Li. 

Comparison of the experimental spectrum of ^B with 
the predictions of ah initio models is shown in Figure [S] 
Clearly, there is no unified picture for the level structure 
of ^B - ®Li isotopes from the available array of ah initio 
calculations. There is a certain general trend, however. 
All of them produce a 0+ as the second or third excited 
state, always below the known 1^ state (experimentally 
found at 3.2 MeV in ^Li) and the 2+ state, found in 
this work at 2.55 MeV. Our experimental result seems to 
confirm this prediction. The 2+ state is generally found 
at higher excitation energy in ah initio calculations than 
observed in this work for '^B. 

In summary, the excitation function for p+^Bc elas- 
tic and inelastic scattering was measured in the energy 
range of 1.6 - 2.8 MeV and the cm. angular range of 
132 - 148 degrees. The R-matrix analysis of the excita- 
tion functions indicates that new low-lying states in ^B 
have to be introduced to explain the large inelastic cross 
section with a well defined peak at 2.5 MeV. These new 
states are suggested to be the 0+ at 1.9 MeV and 2+ at 
2.55 MeV with width 610 keV and 340 keV respectively. 
Evidence for the 2+ state at 2.55 MeV is reliable. The 0+ 
at 1.9 MeV can be considered as tentative due to uncer- 
tainties associated with the coupled-channels calculations 
and possible contributions from the tails of the higher en- 
ergy resonances into the inelastic scattering cross section. 
Accurate measurement of the p-|-^Be excitation function 
for inelastic scattering in the energy range from 0.7 to 
2.0 MeV and in a broad angular range should provide a 
definitive answer on the existence of the 0^ resonance. 

The TDCSM analysis of the ^Be-|-p scattering has been 
performed. The TDCSM reduces the number of free pa- 
rameters in the fit and links directly nuclear structure to 
nuclear reaction cross sections while treating the contin- 
uum self-consistently. The role of TDCSM in our work 
was to provide important guidelines for constraining the 
R-matrix fit. We believe that it is an essential tool for 
analysis of resonance scattering data. 

Comparing new experimental results to the predictions 
of the ah initio models 0, y, [la, [i3| we notice that there 
is no unified picture for the level structure of *B - ^Li 
isotopes. We hope that new experimental data on the 
structure of exotic nuclei (including those presented here) 
will serve as a guide for construction of more accurate ab 
initio models. 
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